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ABSTRACT 



Aims. We search for penumbral magnetic fields of opposite polarity and for their correspondence with downflows. 

Methods. We used spectropolarimetric HINODE data of a spot very close to disk center to suppress the horizontal velocity components 

as much as possible. We focus our study on 3-lobe Stokes V profiles. 

Results. From forward modeling and inversions, we show that 3-lobe profiles testify to the presence of opposite magnetic fields. They 
occur predominately in the mid and outer penumbra and are associated with downflows in the deep layers of the photosphere. 
Conclusions. Standard magnetograms show that only 4% of the penumbral area harbors magnetic fields of opposite polarity. If 3-lobe 
profiles are included in the analysis, this number increases to 17%. 

Key words. Sunspots - Sun: photosphere 



1. Introduction 

Our knowledge of the penumbral magnetic and velocity field 
has increased tremendously (see e.g. Rempel & Schlichenmaier 
2011; Borrero & Ichimoto 2011). The discrepancy between the 
inclination of the average magnetic field and the observed veloc- 
ity field is one example. While a zenith angle of less than 90° 
is deduced for the azimuthally averaged magnetic field (Lites 
et al. 1993), the Evershed flow (Evershed 1909) appears paral- 
lel to the surface throughout the penumbra (Maltby 1964). High 
resolution observation show that the inclination of the penum- 
bral magnetic field is larger within the dark penumbral filaments, 
e.g. Beckers & Schroter (1969); Schmidt et al. (1992); Title et al. 
(1993). Today it is accepted that the inclination of the penumbral 
magnetic field has strong fluctuations on scales of less that 1" 
(Solanki & Montavon 1993; Martinez Pillet 2000; Bellot Rubio 
et al. 2004). It was further argued that the Evershed flow occurs 
along the dark penumbral structures, where the field is more par- 
allel to the solar surface and possibly weaker (Stanchfield et al. 
1997; Schmidt & Schlichenmaier 2000; Langhans et al. 2005). 

However, the correct mode of penumbral magnetoconvection 
is still a controversy. Evidence exists that the upflows in the in- 
ner penumbra are the source of this flow, while the downflows 
in the outer penumbra are the regions where it returns below 
the solar surface (Schlichenmaier & Schmidt 2000; Ichimoto 
et al. 2007; Franz & Schlichenmaier 2009). This is in accordance 
with the 'flux-tube' model (Meyer & Schmidt 1968; Montesinos 
& Thomas 1997; Schlichenmaier et al. 1998) where penumbral 
plasma is channeled along bundles of strong magnetic field. This 
model, in turn, requires a large amount of magnetic field sub- 
merging with these downflows. This should show up as a sig- 
nificant amount of opposite polarity signal in penumbral mag- 
netograms. There have been reports of magnetic fields return- 
ing within the penumbra (Westendorp Plaza et al. 1997, 2001; 
del Toro Iniesta et al. 2001), but even in magnetograms of high 
spatial resolution; e.g., (Langhans et al. 2005) the amount of re- 



turn flux is far less than what is expected from penumbral 'flux 
tube' models (Solanki 2003; Schlichenmaier & Solanki 2003). 
In the debate about the correct mode of magneto convection and 
the penumbral structure, this has been taken as an argument to 
support the 'gappy' penumbral model (Scharmer & Spruit 2006; 
Spruit & Scharmer 2006). However, Ichimoto et al. (2007) used 
HINODE data to compute magnetograms in the far line wing 
and found several opposite polarity patches that are associated 
with penumbral downflows. Nevertheless, detailed studies of the 
amount of penumbral return flux are still missing. 

In this contribution we study 3-lobe Stokes V profiles, which 
indicate the presence of magnetic fields of opposite polarity, if 
the spot is observed at disk center. Considering these profiles, 
we find more opposite polarity than previously reported. We use 
spots that are as close as possible to disk center, and compute 
maps of various physical quantities in Sec. 2. We discuss 3- 
lobe V-profiles that have an extra redshifted lobe in the red wing 
and V-profiles that have a blue-shifted hump in the blue wing. 
We use geometrical models (Sec. 3) and inversion techniques 
(Sec. 4) to model these profiles. We conclude that blue-hump 
profiles are location of upflows with the spot polarity, and 3-lobe 
profiles are locations of opposite polarity. We find that 3-lobe 
profiles preferentially occur at locations of downflows (Sec. 5). 
In Sec. 6 we demonstrate that 'classical' magnetograms detect 
only a small fraction of opposite polarity, since they do not con- 
sider 3-lobe profiles. We quantify how the 3-lobe profiles are 
distributed throughout the penumbra. In Sec. 7 we discuss our 
results. Summary and conclusion are given in Sec. 8. 



2. Observation & data processing 

We used data obtained by the spectropolarimeter (SP) (Lites 
et al. 2001) of the solar optical telescope (SOT) (Tsuneta et al. 
2008) onboard the Hinode satellite. The SP records the Stokes 
spectra of the two iron lines at 630.15 nm and 630.25 nm, with 
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Lande factors of g e ff = 1 .67 and g e fj =2.5, respectively. By scan- 
ning the spectrograph slit across the target in steps of 0."15 and 
with a wavelength sampling of 2.15 pm/pixel, two-dimensional 
maps of the field of view (FOV) were obtained. Since the width 
of the slit is equivalent to 0."16 and the pixel size along the slit 
is 0."16 on average, normal SP scans provide a spatial resolution 
of 0."32 (Centeno et al. 2009). For an exposure time of 4.8 s per 
slit position, the 1 cr noise level in the Stokes spectra is around 
10 3 -I c , where I c is the continuum intensity. 

Table 1. Data samples used in this study. 



Name 


NOAA 


Date of 


Heliocentric 




Active region 


observation 


angle @ [°] 


Spot A 


10923 


Nov 14 lh 2006 


10.8-5.4 


SpotB 


10923 


Nov 14 th 2006 


12.1 - 3.9 


SpotC 


10930 


Dec 1 1 th 2006 


9.8-3.6 


SpotD 


10933 


Jan 05 th 2007 


4.8- 1.5 



For our study, we used the raw data summarized in Table 1 
and reduced it, using the IDL routine 'sp_prep.pro' provided 
with the 'Solar Soft' package. In Sec. 5 and 6 we focus on dataset 
Spot D, because it is the one located closest to disk center. From 
the Stokes spectra, we computed maps of: 

- continuum intensity: I = I c /I qs where I qs is the average con- 
tinuum intensity of the quiet Sun 

- linear polarization: Q(A) = Q(A)/I C and U(A) = U(A)/I C with 
Q(A) and U(A) as the measured Stokes Q and U profiles 

- circular polarization: V(A) = V(A)/I C where V(A) represents 
the measured Stokes V profile 

- total polarization: P tot = / [(Q(A) 2 + U(A) 2 + V(A) 2 )] 1/2 dA 

- net circular polarization 1 : N — JV(A)dA 

- area asymmetry : J{ = N/J |V(A)|dA 

- Doppler velocity along the line of sight: Vdop 

- and a mask of 3-lobe Stokes V profiles: V3_i b e 

The calculation of the maps of I, V, P to t, N and S\ is straight- 
forward. The map of Vd op is constructed from the shift of the line 
wing with respect to a reference wavelength. We refer the reader 
to Franz & Schlichenmaier (2009) for a detailed explanation on 
the calculation of Vd op including an absolute calibration of the 
wavelength scale. 

The mask of 3-lobe Stokes V profiles (V3_i b e ) is constructed 
from the Fe 630.25 nm line, because it is a simple line triplet 
and has the higher magnetic sensitivity. This makes it easier to 
identify the third lobe on the red side of the regular Stokes V 
profile. As indicated by the colored regions in the upper right 
hand panel of Fig. 1, we split the spectrum between 630.205 nm 
and 630.325 nm into a sequence of three windows. The first two 
windows (blue and orange) have a bandwidth of 20 pm, while 
the third one (red) is 45 pm wide. This sequence of spectral win- 
dows is scanned as a whole across the Fe I 630.25 nm line. To 
identify 3-lobe Stokes V profiles, a number of conditions have 
to be fulfilled, a) The lobes in the blue (left) and red (right) win- 
dows are positive, b) The maxima of the lobes and the spectral 
points next to it are positive and above the 3tr noise level, c) The 
lobe of the V profile in the orange (middle) window is negative. 



1 For technical reasons the integral is taken over the entire spectral 
region of the HINODE SP covering both Fe lines. This quantity has 
also been termed broad band polarization (e.g. Illing et al. 1975). 



d) The minimum itself and the spectral measurements next to it 
are negative and above the 3<x noise level. The requirement that 
the neighboring spectral measurements have to be above the 3<x 
threshold and have the same sign as the extrema was applied to 
exclude artifacts caused by the magneto-optical effect, i.e. addi- 
tional peaks around the line core position in Stokes V. A similar 
method for automatic detection of blue-hump profiles encoun- 
ters difficulties and has not been devised. 



3. Blue-hump and 3-lobe Stokes V profiles 

There has been a long history of studies of the asymmetry in 
Stokes V profiles. Babcock (1951) already noted differences in 
intensity and width between the two components of circular po- 
larized light from a variable star. He termed this phenomenon 
the crossover effect 2 . Grigorjev & Katz (1972) were among the 
first to observe Stokes V profiles with more than two lobes in 
the spectra of sunspots. By systematically analyzing the asym- 
metry of Stokes profiles Sanchez Almeida & Lites (1992) and 
Schlichenmaier & Collados (2002) found that Stokes V profiles 
with more than two lobes predominate, but are not exclusive 
(Bellot Rubio et al. 2002), along the magnetic neutral line. 

With HINODE these peculiar profiles are regularly mea- 
sured throughout the entire penumbra. Ichimoto et al. (2007) 
and Franz & Schlichenmaier (2010) found that Stokes V profiles 
from penumbral upflow regions contain an additional hump on 
the blue side of an antisymmetric Stokes V profile (blue-hump 
profiles), while the corresponding profiles from downflow re- 
gions contain a third lobe of opposite sign on the red side (3-lobe 
profiles) - cf. Fig. 1. 
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Fig. 1. Stokes V profiles from typical penumbral up- and down- 
flows observed with HINODE (crosses). Top row: Synthetic pro- 
files (solid) radiative transfer calculations within a two-layer at- 
mosphere. The colored regions in the top right plot indicate the 
spectral regions that are used to identify 3-lobe profiles. Bottom 
row: Synthetic spectra (solid) originating in a two-component at- 
mosphere. 

There are two simple atmospheric configurations that yield 
blue -hump or 3-lobe Stokes V profiles. The first is a two-layer 
atmosphere. It includes strong jumps (or gradients) of the at- 
mospheric parameters along the line of sight (LOS) and within 
the line-forming region, e.g. Auer & Heasley (1978) or Landi 



2 We use the term 3-lobe profiles instead of crossover effect to differ- 
entiate these profiles from stellar profiles. 
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Degl'Innocenti & Landolfi (1983). The second configuration is 
a two-component atmosphere. It mixes two laterally displaced 
but spatially unresolved atmospheres (in these atmospheres no 
gradient of Doppler velocity, strength, inclination, and azimuth 
of the magnetic field are present along the LOS) within one res- 
olution element, e.g. Bellot Rubio et al. (2004). 

3. 1 . Geometrical models 

We studied the ability of these two configurations to repro- 
duce blue-hump or 3-lobe Stokes V profiles by means of ra- 
diative transfer calculations. We used the synthesis module of 
SIR (Stokes Inversion based on Response functions Ruiz Cobo 
& del Toro Iniesta 1992), together with a modified version 
of the Harvard Smithsonian Reference Atmosphere (HSRA) of 
Gingerich et al. (1971). The changes we introduced in the HSRA 
to investigate the two-layer atmosphere are listed in Table 2. 



Table 2. Two-layer atmosphere: Modification to HSRA. 



Optical depth 


Field strength Zenith angle 


Velocity 


log(T 50 o) 


B [G] y [ ] 


Vdop [km s -1 ] 


Upflow case 


0.0 to -0.5 


1000 60 


-6.5 


-0.6 to -3.0 


1500 20 


0.0 


Downflow case 


0.0 to -0.5 


1000 120 


8.5 


-0.6 to -3.0 


700 60 


0.0 



To mimic a two-component atmosphere, we performed radia- 
tive transfer calculations within a modified version (cf. Table 3.) 
of the HSRA. To obtain the final Stokes V profile, we mixed the 
two Stokes V profiles from the individual atmospheric compo- 
nents with different filling factors. 



Table 3. Two-component atmosphere: Modification to HSRA. 



Filling factor of 


Field strength Zenith angle 


Velocity 


component [%] 


B [G] y [°] 


v dop [km s" 1 ] 


Upflow case 


30 


1000 60 


-6.5 


70 


1500 20 


0.0 


Downflow case 


30 


1000 120 


8.5 


70 


700 60 


0.0 



Figure 1 summarizes the outcome of our experiment. It 
can be seen that both the two-layer atmosphere and the two- 
component atmosphere is able to reproduce blue-hump and 3- 
lobe Stokes V profiles. Keep in mind that these models are con- 
structed to reproduce distinct features and not to minimize the 
difference between the observed and synthetic profiles. The lat- 
ter is done in Sec. 4 by applying an inversion method. 

We found a vast range of possible atmospheric configura- 
tions that yield blue-hump or 3-lobe Stokes V profiles. However, 
in all these realizations it is necessary to introduce a) magnetic 
fields of sufficient strength in both components or layers, b) an 
unshifted and a blueshifted component or layer for the upflow 
case, c) an unshifted and a redshifted component or layer for the 



downflow case, d) magnetic fields of the same polarity in the 
blueshifted component or layer, and e) magnetic fields of oppo- 
site polarity in the redshifted component or layer. 

In both configurations magnetic fields of opposite polarity 
are present in downflow regions. In the two-component atmo- 
sphere they simply lie next to each other in one resolution el- 
ement. In the two-layer atmosphere the downflows are concen- 
trated in the lower layers of the photosphere, < log(r) < -0.5, 
where the magnetic field returns below the solar surface, y > 
90°. In the layers above, 0.6 < log(t) < -3.0, the atmosphere is 
at rest and the magnetic field has the same polarity as the umbra. 

3.2. Asymmetry of blue-hump and 3-lobe profiles 

So far our results have not allowed determining whether gra- 
dients along the LOS or laterally unresolved structures within a 
resolution element are at the origin of blue-hump and 3-lobe pro- 
files. However, it is possible to use the net circular polarization 
(AO or the area asymmetry 3 (^) to distinguish between these 
two possibilities. According to Auer & Heasley (1978), a gra- 
dient of Vdop along the LOS is not only necessary, but sufficient 
to produce non-vanishing values of N, hence 3L If a gradient in 
Vdop is present, gradients of other parameters, i.e. the strength, 
the inclination or the azimuth of the magnetic field, alter the 
amount of J\ (Landolfi & Landi Degl'Innocenti 1996). 

For the two profiles presented before we measure ^°p S = 

19% and &f^ n = 84%, respectively. Such high values of J[ are 
common for Stokes V profiles from penumbral up- or downflows 
and need to be reproduced by the synthetic profiles. 

The two-layer atmosphere produces a significant amount of 
area asymmetry: jTlup 61 = -24% and ^down = 24%. However, 
the synthetic upflow profile yields the wrong sign of J{ when 
compared to observations. 

The two-component atmosphere is unable to produce any 
area asymmetry, since the atmospheric parameters do not change 
along the LOS. Therefore, we do not consider this configuration 
in the following. 

4. Inversion 

We used the inversion tool of the SIR code (Ruiz Cobo & del 
Toro Iniesta 1992) to create an atmosphere as simple as possi- 
ble, but still capable of reproduceing all spectral features (blue- 
hump and 3-lobe), as well as the area asymmetry of the measured 
Stokes V profiles. In contrast to forward modeling, this tech- 
nique minimizes the differences between observed and synthetic 
Stokes profiles, yielding a maximal goodness of the fit. In SIR, 
physical parameters are interpolated between 'nodes', which are 
placed in the atmosphere. In the iteration process, the values at 
the nodes are changed. If there is only one node for a parameter, 
its initial function along the line-of-sight is shifted according to 
the change in the node. When the iteration process converges, a 
subsequent iteration loop (i.e. cycle), with an increased number 
of nodes, can be processed. 

4.1. Setup and procedure 

We performed a two-cycle inversion and used a modified version 
of the HSRA as an initial model. We randomized the temperature 

3 In contrast to N, is normalized to the area of the profile. This 
makes the latter more appropriate when comparing the asymmetry of 
Stokes V profiles with different amplitudes. 
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Fig. 2. Top part clockwise: Representative Stokes I, Q, V, and 
U profiles (black) from an upflow region. The red lines indicate 
the result of radiative transfer calculations within the atmosphere 
plotted below and the dashed vertical lines represent the respec- 
tive vacuum wavelength of the transition. Bottom part clock- 
wise: Temperature, Doppler velocity, zenith angle, and magnetic 
field strength (black) with errors (shaded gray) at different opti- 
cal depth values. The dashed lines indicate the temperature strat- 
ification in the HSRA model, plasma at rest, and a zenith angle 
of 90° in the respective panels. 

(T) distribution around the values of that of the HSRA, while 
the values for the electron pressure (p e ) and microturbolence 
(Vmic) were kept. We added random values of Doppler velocity 
(vdopX magnetic field strength (B), zenith angle (y), and azimuth 
((b), but ones that are constant in optical depth (t). To ensure 
that the final atmosphere does not represent a local minimum 
of the x 2 manifold 4 , the inversion of a given Stokes vector was 
performed 100 times with randomized input model atmospheres. 
From the resulting pool of solutions, the fit with a minimal x 2 
was selected. 

In the first cycle, we used one node in all atmospheric pa- 
rameters except for p e -, which remains the same as in HSRA. 
The macroturbolence (v mac ) is also a free parameter, but we did 
not consider a straylight component. This leads to seven degrees 
of freedom (DOF). 

In the second cycle, we allowed for three nodes in the strat- 
ification of T and two in that of B and y, while the rest was not 
changed. For Vd op , it was necessary to allow for seven nodes to 
be able to retrieve the shape of the observed Stokes V profiles. 
In fact, it is impossible to model the additional hump in upflow 
profiles with less DOF in Vd op - This is feasible for the red-shifted 

4 X 2 is a weighted difference between all observed and synthetic 
Stokes parameters divided by the number of free parameters. 



Fig. 3. Same as Fig. 2 but for a representative profiles from a 
downflow region. The magnetic field changes its polarity in the 
deep photospheric layers. 



profiles, but only at the cost of increasing the number of nodes 
in B and y to three. The reason for this lies in the mode of op- 
eration of the code itself. SIR distributes the nodes equally with 
t, which makes it difficult to model a flow mainly present in the 
deep atmospheric layers. All together there are 17 DOF in the 
second cycle. 

4.2. Upflow atmosphere 

Figure 2 represents the solution for one upflow case. The four 
panels at the bottom illustrate the stratification of the atmo- 
spheric parameters. Depicted are T, Vd op , B, and y mag only, be- 
cause the other parameters where either not inverted (p e ) or re- 
main constant throughout the line-forming region (v m [ c , v mac , 
and dp). The variation in the atmospheric parameters along the 
LOS is plotted between log(r) = -3 and log(T) = 0, where the 
sensitivity of the two lines is highest. 

In this model, T increases monotonously from 4300 K at 
log(t) = -3 to 6450 K at log(T) = 0, with an error of +30 K. 
Except for a little plateau around log(x) = -2, this distribution 
is similar to that of the HSRA (cf. dashed line in the respec- 
tive plot of Fig. 2). The strongest flows occur in the low lay- 
ers of the atmosphere, where the largest amplitudes of Vd op are 
present. Between log(r) = -0.7 and the surface, the amplitude 
increases from Vd op = km s _1 to Vd op = -7 km s while, 
within the uncertainties, almost no Vd op is measured in the layers 
above (-3.0 < log(x) < -0.7). B and y mag remain almost con- 
stant throughout the atmosphere with values of 60° +5.5° and 
1 .5 ± 0.16 kG respectively. 
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Fig. 4. Top: Average Stokes V profiles of different bins of Doppler velocity for the penumbra of dataset Spot D. White represents 
the positive lobe, black the negative lobe. The picture is saturated at a continuum intensity of 2.5% and all measurements below the 
3-cr noise level were artificially set to zero. Bottom: Examples of averaged profiles in various velocity bins. 



4.3. Downflow atmosphere 

The solution for one downflow case is depicted in Fig. 3. When 
comparing these profiles to the upflow case (cf. Fig. 2), it be- 
comes apparent that Stokes I is shifted towards longer wave- 
lengths, especially in the wing of the line. The value of I c is 
below I qs , and the maximal amplitudes of Stokes Q, U, and V 
are lower than in the upflow profile. 

In the resulting model, T increases monotonously from 
4350 K at log(x) = -3 to 5300 K at log(T) = 0, with an er- 
ror of roughly +40 K. At lower layers, this atmosphere is cooler 
than the HSRA (Ti og(T) (-3) = 4400 K and T log(T) (0) = 6400 K), 
which explains the 25% drop in I c . The distribution of T shows 
a little plateau around log(x) = -1, a feature not present in 
the HSRA. For Vd op , strong flows are found in the low lay- 
ers of the atmosphere. Their amplitude increases monotonously 
from Vdop = km s at log(T) = -1.5 to high values of 
Vdop = 8.5 km s at the surface. There is also a slight downflow 
in the higher layers (Vd op = 1-5 km s _1 for -3 < log(T) < -2.5), 
but no flow for -2.5 < log(T) < -1.5. The uncertainty in Vd op is 
about +0.5 km s _1 . B increases from 1.1 kG to 1.6 kG between 
log(T) = -3 and log(T) = 0. The relatively large error of +0.35 
kG makes it impossible to rule out a constant B throughout the 
atmosphere. In contrast to the upflow example, y increases from 
25° at log(T) = -3 to 130° at log(T) = 0, with an error of ±10°. 
This means that the magnetic field reverses its polarity some- 
where between log(x) = -1.5 and log(T) = -0.9, an atmospheric 
height reached by the strong downflows of the lower layers. 

4.4. Discussion 

The area asymmetry of the resulting blue-hump profiles from an 
upflow region has an area asymmetry of = 7% (compared 

toy[* s = 19%)andQ up = 0.89. For the 3-lobe profile, we obtain 

^down = 89% (compared to fl*^ = 84%), while Q do „n = 0.78. 
The inverted profiles that reproduce the observed features (blue- 
hump and 3-lobe). Differences exist in the area asymmetry, es- 
pecially for the blue-hump profiles. This is because the SIR code 
minimize the difference between the observed and the synthetic 
Stokes profiles, and not between observed and synthetic M. 



5. From blue-hump to 3-lobe Stokes V profiles 

So far we have focused on one example profile for the up- 
and downflow cases, respectively. In the following we study 
all Stokes V profiles, and we demonstrate that the occurrence 
of 3-lobe profiles increases continuously with Doppler velocity 
(Vdop)- To this end we focus on dataset Spot D, because it is the 
one located closest to disk center and penumbral up- and down- 
flows are most visible. 

5. 1 . Behavior of average profiles of velocity bins 

For this investigation we ordered the penumbral Stokes V pro- 
files according to their Doppler velocity. Average profiles were 
computed for velocity bins of 5 m s _1 , e.g. 750 m s _I < Vd op < 
755 m s _1 . The result is depicted in Fig. 4, where the top 
panel shows the averaged profiles in all velocity bins between 
±1.5 km s _1 . The five plots - a) through e) - in the bottom row 
represent examples of different velocity bins. We set all measure- 
ments below the 3-ff noise level to zero, and omitted all velocity 
bins in which less than five Stokes V profiles were averaged. 

Figure 4 shows that the white area (the positive lobe) is wider 
than the black area (the negative lobe) in velocity bins represent- 
ing strong upflows. This is due to the additional hump on the blue 
side of the regular Stokes V profile (cf. Plot a) in Fig. 4. With de- 
creasing upflow velocity, the width of the positive lobe decreases 
continuously as the additional hump disappears (Fig. 4 b). At 
zero Doppler shift, both lobes have the same width, and the pro- 
file appears antisymmetric (Fig. 4 c). For velocity bins represent- 
ing downflows greater than 0.6 km s we detect an additional 
lobe above the 3-cr noise level in the averaged profile around 
+25 pm from the linecore of the Fe I 630.25 nm line. With in- 
creasing downflow velocity, this additional lobe grows in width 
and amplitude (Fig. 4 e). 

Since the penumbra in Spot A, B, and C is located farther 
away from disk center, the projection of the horizontal Evershed 
outflows onto the line of sight yields Doppler shifts as well. 
Therefore, we conducted our analysis only in an area defined 
by a cone angle of ±30° perpendicular to the line of symmetry - 
i.e. a line connecting the center of the spot with the center of the 
solar disk. In these regions the contribution of the radial outflows 
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Fig. 5. Blue crosses and left axis: Percentage of 3-lobe profiles in 
bins of increasing downflow velocity for Spot D. Red histogram 
and right axis: Fraction of profiles in velocity bin compared to 
the overall number of penumbral profiles. Solid blue and solid 
red lines: A second order polynomial and an exponential decay 
fit, respectively. 



to the velocity signal is always less than 50% when compared to 
regions along the line of symmetry. If we correct further for the 
different spot polarities in the respective datasets, Spot A, B, and 
C (not shown here) yield the same pattern of averaged Stokes V 
profiles as presented in the top panel of Fig. 4. 

5.2. Fraction of 3-lobe profiles in velocity bins 

Finally we want to exclude the possibility that the additional lobe 
in the profiles - panels a) through e) in Fig. 4 - is an artifact 
of the averaging process. We show that the additional lobe in 
the depicted profiles is not caused by a superposition of regular 
Stokes V profiles with various line shifts and polarities, but are 
due to an actual increase in the number of 3-lobe profiles in the 
respective velocity bin. 

We use the profiles of the downflow regions, because it is 
easier to detect the additional lobe above the 3-cr noise level on 
the red side of the profile. In each velocity bin we counted the 
number of profiles that obey the criteria for 3-lobe profiles - 
see Sec. 2 - and computed their percentage on the total number 
of profiles. We omitted all bins with fewer than five profiles. 

Figure 5 demonstrates that the number of profiles in the re- 
spective bin decreases exponentially with velocity (cf. Fig. 3 in 
Franz & Schlichenmaier 2009), while the percentage of 3-lobe 
profiles on the overall number of Stokes V profiles in the respec- 
tive velocity bin increases. 

For 0.0 km s _1 < v dop < 0.2 km s _1 10% to 20% of all pro- 
files in a velocity are 3-lobe profiles. This percentage increases 
almost linearly until Vd op = 0.5 km s and already for Vd op = 
0.55 km s _1 , individual velocity bins exist in which more than 
50% of all profiles contain a third lobe. For higher downflow ve- 
locity, the scatter increases significantly, primarily because the 
number of profiles in a single bin decreases. According to the 
polynomial fit, 3-lobe profiles dominate for Vd op > 0.7 km s _1 
and saturate around 70% for Vd op > 1 .0 km s . 

Thus, it is indeed the number of 3-lobe profiles in one ve- 
locity bin that increases with downflow velocity and causes 
the additional lobe in the average profiles of Fig. 4. Note that 
3-lobe profiles appear in the bins of small downflow velocity 



(0.0 km s _1 < Vdop ^ 0.5 km s ) and sometimes even in bins 
of upflows. This could be due to the spatial resolution of the in- 
strument, which can result in a mixing of blue-hump and 3-lobe 
profiles from nearby up- and downflows, cf. Sec. 7. 



6. Distribution of penumbral 3-lobe profiles 

Magnetic fields of opposite polarity are commonly inferred from 
'normal' 2-lobe Stokes V profiles of opposite sign. Here, we 
demonstrated that opposite polarity in a particular pixel can also 
be inferred from the presence of 3-lobe profiles (if the spots is at 
disk center), leading to a much larger fraction of opposite polar- 
ities in a penumbra. This is of special relevance for penumbral 
models: Spruit & Scharmer (2006), for example, argue against 
'flux tube' models stating that these models need a large amount 
of submerging magnetic flux, which is not (or only partially) 
detected in magnetograms by Langhans et al. (2005) (cf. also 
Solanki & Schmidt 1993; Westendorp Plaza et al. 1997). 

Here we show that there is much more return flux than de- 
tected with 'Langhans' -magnetograms, if one takes into account 
3-lobe profiles. In Sec. 6.1 we use our data to construct (i) mag- 
netograms mimicking Langhans - cf. 6.1, (ii) compute magne- 
togram masks of V-profiles with opposite sign following Lites 
- cf. Sec. 6.2, (iii) and merge the latter mask with a mask of 
3-lobe profiles to quantify all locations of opposite polarity - 
cf. Sect. 6.3. Finally (Sec. 6.4) we compare our mask to a map, 
in which only the V-signal in the far red wing is displayed. 

6. 1 . Simulated SOUP@SST magnetogram 

The study of Langhans et al. (2005) was conducted at the 
Swedish Solar Tower (SST), using the Lockheed Solar Optical 
Universal Polarimeter (SOUP). This Lyot-type filtergraph mea- 
sures the intensities of left and right circular polarized light with 
a bandpass of FWHM 7.2 pm at various wavelengths across 
the spectral line (Title & Rosenberg 1981). In the case of Fe I 
630.25, 'Langhans' -magnetograms are constructed like 

jLCP _ jRCP 



with fr. and Ifr being the intensities of left and right circular 

blue blue & & 

polarized light, measured by SOUP detuned -5 pm to the blue 
with respect to the linecore. 

To simulate SST observations with the SOUP instrument, we 
first convolved the respective HINODE Stokes V measurements 
with a Gauss profile (100% transmission at -5 pm from the line 
core of Fe I 630.25 nm with a FWHM of 7.2 pm). In a sec- 
ond step we integrated the obtained profiles over the full spectral 
range of SOUP and assigned a magnetic field of opposite polar- 
ity to the respective pixel if the result is the opposite sign to the 
umbral pixels. 

Masking opposite polarity pixels with a red color, the mor- 
phology of this map is displayed in the upper left hand panel 
of Fig. 6. We see that opposite polarity pixels are predomi- 
nantly detected close to the outer penumbral boundary. Only 
2% (1506 out of 74842 pixels) of the penumbral area is of op- 
posite polarity. From these, 91% are locations of downflows 
(Vdop > 0.1 km s _1 ). From all downflow pixels, 5% (1350 out 
of 26064 pixels) have opposite polarity. 
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Fig. 6. Clockwise from top left: Penumbral magnetic fields of opposite polarity as they appear in a simulated SOUP@SST magne- 
togram, in a standard HINODE SP magnetogram, in a HINODE magnetogram computed in the far line wing of Fe I 630.25 nm, 
and in a standard HINODE magnetogram if 3-lobe Stokes V profiles are taken into account as well. 



6.2. Standard SP magnetogram 



Using the standard SP data pipeline (SSWIDL), maps of the 
apparent longitudinal magnetic field, B??, can be obtained. 



According to Lites et al. (1999, 2008) 5 - the apparent longitudi- 
nal magnetic field is proportional to 



gapp 
Ion 



/ A ;V(A)dA-/ A ;V(A)dA 
Ic£dA 



(2) 



5 Even though the definition from Bj? in these publications is 
slightly different to Eq. 2, they yield identical maps of B?? in the 
penumbra. Differences occur in the signal-to-noise ratio of weak Stokes 
V signals, which are not very common in the penumbra. 
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where I c denotes the continuum intensity and Ao the zero- 
crossing wavelength. The upper and lower limits of integration 
are represented by A r and Ab, which are set to A r /b = Ao+30 pm, 
respectively. The polarity of the flux in the maps of Bj^ is thus 
determined from the sign of the area of the two lobes around the 
zero crossing. We manipulated the IDL routines in such a way 
that only the Fe I 630.25 nm line is considered. 

The corresponding 'SP standard' -mask of opposite polar- 
ity is shown in the upper right hand panel of Fig. 6. As in the 
'Langhans'-mask, the pixels are located in the outer penumbra, 
but almost two times as many are found with this method: 4% 
(2783 out 74842 pixels). From these, 92% correspond to down- 
flows. All opposite polarity pixels of the 'Langhans'-mask are 
contained in this 'SP standard' -mask. Now 10% (2539 out of 
26064 pixels) of all downflow pixels have opposite polarities. 

6.3. Standard SP magnetogram and 3-lobe profiles 

In this study we determine the amount of opposite polarities in 
the penumbra including 3-lobe Stokes V profiles. To this end we 
merge the 'SP standard' mask with the mask of 3-lobe profiles 
(V3_i be)- This is justified since regular but mirrored Stokes V 
profiles remain undetected in the V3_i he mask (cf. the criteria 
defined in Sec. 2), but they are still tracers of opposite polarity. 

Interestingly, the two sets are not disjunct. We find that only 
half of the profiles in the 'SP standard' mask are mirrored Stokes 
V profiles. The other half contains 3-lobe profiles in which the 
additional lobe on the red side becomes larger than the difference 
in the areas of the two regular lobes. If Eq. 2 is applied to the 
latter group of profiles, the result will be interpreted as magnetic 
fields of opposite polarity. 

The superposition of V3_i b e and the 'SP standard' mask is 
shown in the bottom left hand panel of Fig. 6. Now 17% (12594 
out of 74842 pixels) of the penumbral area shows magnetic fields 
of opposite polarity and 87% of these opposite polarities are 
located in downflows. Thus, the amount of all the penumbral 
downflows that harbor a magnetic field of opposite polarity in- 
creases to 40% (10264 out of 26064 pixels). 

In contrast to the two previous masks, the change in the mor- 
phology is striking - cf. lower left hand panel in Fig. 6. The 
return flux is very prominent in the mid and outer penumbra, but 
occurs in the inner penumbra too. This is also the case in the 
penumbral side pointing towards disk center. This implies that 
the 3-lobe profiles we use as tracers of penumbral return flux 
cannot be due solely to horizontal Evershed motions. 

6.4. Magnetogram in the far line wing 

Ichimoto et al. (2007) computed magnetograms in the far red 
wing of the Fe 630.15 nm line to illustrate penumbral magnetic 
fields of opposite polarity. We followed this procedure, and com- 
puted a magnetogram at +20.4 pm from the line core of the Fe 
630.25 nm line in order to compare such a measurement to our 
analysis. In contrast to Ichimoto et al. (2007), we considered a 
Stokes V profiles as an indicator of opposite polarity only if the 
signal at +20.4 ± 2.15 pm, i.e the neighboring spectral measure- 
ments, are above the 3 cr noise level, too. 

The lower right hand panel of Fig. 6 shows a mask of oppo- 
site polarities in the far wing magnetogram. They occupy 13% 
(9806 out of 74842 pixels) of the penumbral area and 85% of 
them are located in downflows. As a result the fraction of penum- 
bral downflows with magnetic fields of opposite polarity is 32% 
(8364 out of 26064 pixels). Overall, the opposite polarity patches 



are smaller when compared to the previous mask, i.e. left hand 
panel of Fig. 6, but the overall morphology is very similar. 

Our far wing magnetogram recovers 86% of the opposite po- 
larity pixels of the 'SP standard' -magnetogram. Differences oc- 
cur for example if the third lobe of a 3-lobe profile is not above 
the noise level at the specified wavelength. Furthermore, there 
are Stokes V profiles with only one or sometimes even with four 
lobes. We cannot rule out that such profiles yield a opposite po- 
larity signal in the 'SP standard' magnetogram, but at the same 
time, are not detected in the far wing magnetogram. 

Table 4. Penumbral return flux in different magnetograms. 



Method of Penumbral area with Fraction of return flux 

detection opposite polarity[%] in downflows 6 [%] 

Intensity at ^ 5 

fixed wavelength 

Standard SP „ , „ 

4 10 

magnetogram 

Far Wing ^ 
magnetogram 

SP Magnetogram „ ^ 

and 3-lobe profiles 



7. Discussion 

The existence of penumbral downflows harboring magnetic 
fields of opposite polarity is an essential finding in the context 
of penumbral models. That these downflows, together with the 
opposite polarities, occur mainly in the middle and outer penum- 
bra and that (at HINODE resolution) they appear patchy instead 
of elongated, supports the idea of flows in arching magnetic flux 
tubes. 

Such flux tubes are a central ingredient of flux tube models 
(cf. e.g. Meyer & Schmidt 1968; Montesinos & Thomas 1997; 
Schlichenmaier 2002). In the framework of the gappy penumbral 
model (Scharmer & Spruit 2006; Spruit & Scharmer 2006), the 
existence of opposite polarities as well their observed morphol- 
ogy cannot be explained. 

Rempel (2012) argues that the majority of penumbral down- 
flows advect magnetic field, which are still connected to the up- 
per boundary of the simulation box, below the surface without 
changing its polarity. In any case, our findings prove that a large 
fraction of the emerging magnetic flux submerges within the 
penumbra. 

Return flux: We speculate that blue-hump and 3-lobe profiles 
indicate the presence of low-lying loops in the penumbra: the 
upflow footpoint is seen as a blue-hump profile, while the down- 
flow footpoint (of opposite polarity) is seen as a 3-lobe profile. 
This is sketched in Fig. 7. In this scenario, out of the 17% of op- 
posite polarity, 15% exist in the deep line-forming region only 
(3-lobe profiles), while 2% fill the entire line-forming region 
(mirrored profiles). If we assume further that all the submerg- 
ing magnetic flux has a merging counterpart (blue-hump profile), 

6 We consider only penumbral profiles with a downflow signal above 
the detection limit of 0.1 km s _1 , cf. Franz & Schlichenmaier (2009). 
For Spot D these are 26064 out of 74842 penumbral profiles. 
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2-15% = 30% of the penumbral pixels belong to magnetic field 
that returns in the deep photosphere. In this picture 2 -2% = 4% 
of the pixels show arches with apices in higher photospheric lay- 
ers, and 66% of the pixels harbor magnetic fields that continue 
beyond the line-forming region. 




V-profiles normal blue-hump 3-lobe mirrored 



Fig. 7. Proposed configuration of the magnetic field in the line- 
forming layer of Fe I 630.25 nm. 



3-lobe profiles and upflows: Interesting is that 3-lobe profiles 
sometimes exist in upflows (Vd op < -0.1 km s _1 ) as well. These 
profiles look like a superposition of blue-hump and 3-lobe pro- 
files. In our analysis these profiles occupy 1.3% of the penum- 
bral area, and they primarily occur between up- and downflows 
or they surround downflow patches. This suggests that such pro- 
files are a mixture of up- and downflows that are not spatially 
resolved. 

Red-hump profiles: Katsukawa & Jurcak (2010) found an- 
other type of asymmetric Stokes V profiles which show a hump 
on the red side of a regular V profile. These profiles would hint 
at magnetic fields of the same polarity in the downflow region. 
However, patches of red-hump profiles are rather small (w 0."5) 
and far less abundant than the regions of 3-lobe profiles. Since 
these profiles occur co-spatially with a brightening in Ca II H, 
Katsukawa & Jurcak (2010) speculate that such red-hump pro- 
files could be signs of magnetic reconnection. 

Micro Structured Magnetic Atmospheres (MISMA): 

Sanchez Almeida & Ichimoto (2009) studied 3-lobe Stokes V 
profiles under the MISMA hypothesis (cf. also Sanchez Almeida 
2005). Similar to our results, these authors identify penumbral 
magnetic fields of opposite polarity using 3-lobe profiles. They 
find that a multiple component atmosphere is needed to properly 
model these asymmetric profiles. However, unlike our two-layer 
atmosphere, the MISMA hypothesis assumes that the solar at- 
mosphere is structured on scales below the resolution limit of 
any observation available today (Sanchez Almeida et al. 1996). 
Even though our findings are not at odds with the MISMA hy- 
pothesis, we favor an atmospheric model that is as simple as 
possible for explaining all (current) observational features, i.e. a 
two-layer atmosphere. 

Dark filaments and dark cores: Sanchez Almeida & Ichimoto 
(2009) suggest that the weakening of the magnetogram signal 
in dark cores (Bellot Rubio et al. 2007; van Noort & Rouppe 
van der Voort 2008) is due to the mixture of small-scale mag- 
netic fields of opposite polarity within penumbral filaments, 
cf. Sanchez Almeida (2010). We do not find a significant corre- 



lation between opposite polarity (3-lobe Stokes V profiles) and 
dark-cored filaments using HINODE SP data of a penumbra at 
disk center. Only in spots off disk center, e.g. NOAA 10933 ob- 
served on January 7, 2007 at « 20°, such profiles show up in 
dark cores in the limb-side penumbra close to the apparent mag- 
netic neutral line. Such behavior, however, is readily explained 
by the inclination of the flow channel (dark core) under the pro- 
jection of the heliocentric angle. While the observed inclination 
(zenith plus heliocentric angle) of the magnetic field within the 
flow channel is larger than 90°, this is not the case for the back- 
ground atmosphere, which has a zenith angle smaller than 90°. 
Thus, it is a projection effect that yields 3-lobe Stokes V profiles, 
but this does not indicate opposite polarities in the local refer- 
ence frame. Because of the arguments made above, we follow 
the common interpretation that the weaker magnetogram signal 
of the penumbral 'intra-spines' is due to the larger zenith angle 
of the magnetic field, causing a reduction in the Stokes V signal. 

Principal component analysis: Our simple approach of count- 
ing the lobes of Stokes V profiles is limited by the signal-to- 
noise ratio of the HINODE SP. It was shown by del Toro Iniesta 
& Lopez Ariste (2003) that all Stokes profiles can be interpreted 
as an infinite linear combination of a specific set of Hermite 
polynomials which form an orthonormal basis of squared in- 
tegrable functions. Such a basis is noise free, and errors occur 
only because, in practice, an infinite number of terms cannot be 
considered in the series expansion. Using such a principal com- 
ponent analysis, it might be possible to detect an even larger 
number of magnetic fields with opposite polarity. 



8. Summary and conclusion 

We explored strongly asymmetric Stokes V profiles in the 
penumbra of sunspots close to disk center, by using data of high 
spatial and spectral resolution from the Hinode spectropolarime- 
ter. We can confirm earlier findings by Ichimoto et al. (2007) 
that penumbra up- and downflows are accompanied by strongly 
asymmetric Stokes V profiles. The profiles of upflows show an 
additional hump on the blue side, while the profiles from down- 
flow regions contain an additional lobe on the red side, i.e., have 
three lobes. 

Using the 'synthesize module' of SIR, we modeled two sim- 
ple configurations to emulate these asymmetries: Two compo- 
nents that are separated laterally (two-component) or vertically 
(two-layer). Both models are capable of qualitatively reproduce- 
ing the blue-hump and 3-lobe V profiles. For the 3-lobe profiles 
both models need magnetic fields of opposite polarity. However, 
only a vertical separation can produce an area asymmetry in V 
profiles. To quantitatively fit the profiles we performed an in- 
version with SIR. To minimize the degrees of freedom, we used 
a one-component model with gradients along the line-of-sight, 
which provides similar ingredients as the two-layer model to re- 
produce the area asymmetries. We found that the proper repro- 
duction of a 3-lobed profile, and its area asymmetry requires a 
strong magnetic field of opposite polarity and a downflow, both 
concentrated in the deep photosphere (0 < log(x) < —1.5). 

Studying the shape of Stokes V profiles, we ordered them 
in bins of Doppler velocity, and calculated average profiles for 
each bin. From strong blue-shifts to strong red-shifts, the profiles 
transform smoothly from blue-hump profiles to 3-lobe profiles. 
This demonstrates that blue-humps correspond to upflows, while 
3-lobe profiles correspond to s. 
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In 'classical' magnetograms only 4% of the pixels are of op- 
posite polarity; i.e., these magnetograms miss a large fraction 
of opposite polarity. Merging the masks deduced from opposite 
sign of V and 3-lobed shapes, 17% of all penumbral pixels are 
of opposite polarity. From these, 87% correspond to downflows. 

Our numbers are only lower limits. With higher spatial res- 
olution and a better signal-to-noise ratio, we expect that a larger 
amount of penumbral return flux to be detected. 
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